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RECTANGULAR- TUBE SECTIONS - -- 

By Elbridge Z. Stowell and Eugene E. Lundquist 

SUMMARY 



Charts are presented for the coefficients in the for- 
mulas for the critical compressive stress at vrhich cross- 
sectional distortion begins in thin-wall columns of I- , 
Z-, channel, and rectangular-tube sections. The energy "" 
method of Timoshenko was used in the theoretical calcula- 
tions required for the construction of the charts.- The 
deflection equations were carefully selected to give good" 
accuracy. -" ' 

The calculation of the critical compressive stress at 
stresses above the elastic range is briefly discussed. In 
order to demonstrate the use. of the formulas and the charts 
in engineering calculations, two illustrative problems are 
included. 



INTRODUCTION 



In the. design of compression members for aircraft, 
whether'they be stiffeners in stressed-skin structures or 
struts in trussed structures, the allowable stress for the 
member is equal to the. lowest strength corresponding to 
any of the possible types of failure. In references 1 and 
2, all types of column failure are classed under two head- 
ings: 

(a) Primary, or general, failure. 

(b) Secondary, or local, failure. 

Primary, or general, failure of. a column is defined 
as any type of failure in which the cross sections are 
translated, rotated, or both translated and rotated but 
not distorted in their own planes , (fig. l). Secondary, or 
local, failure of a column is defined as any type of fail- 
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ure in which, the cross sections are distorted in their own 
planes hut not translated or rotated (fig. 2). Consider- 
ation is given in this paper- only to local failure. 

One of the factors to he considered in a study of lo- 
cal failure is the dritical compressive stress at which 
the cross section begins to distort. This critical stress 
can usually "be given in coefficient form. In two previous 
papers, coefficients are given for the rectangular tube 
(reference 3) and for columns of channel section and Z- 
section (reference 4). The purpose of the present paper 
is tt> summarize the results embodied in references 3 and 4 
and also to present co'ef f i ci ent a that will permit the crit 
i cal compressive stress to he computed fjor an I-section 
column. 

The calculations required to evaluate the coefficient 
plotted in the charts were made by the energy method of 
Timoshenko (reference 5). Because the calculations are 
long and were made as a part of a more extended study of 
local failure in thin-metal columns, they have been omit- 
.ted from this paper. 

Bernard Rubenstein, formerly of the N.A.C.A. staff, 
performed a large part of the mathematical derivations 
required for the preparation of this paper. 



CHARTS 



The calculation of-the critical compressive stress at 
which cross-sectional distortion begins in any thin-wall 
column is, in reality, a problem in the buckling of thin 
plates, proper consideration being given to the interac- 
tion between the adjacent plates composing the cross sec- 
tion. Each of the sections considered in this paper con- 
sists of two basic plate elements. I'or example, in the— I — 
section these elements are the flange plate and the web 
plate. 

Timoshenko gives the critical stress for a rectangu- 
lar plate under edge compression in the following form 
(reference 6, p.. 605): 
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where 

E is tension-compression modulus of elasticity 
for the material. 

p., Poisson's ratio for the material. 

t, thickness of the plate. 

b, - width of the plate. 

k, a nondimensional coefficient dependent upon 
the conditions of edge support and thedi- 
mensions of the plate. 

This equation can "be used to calculate the critical 
compressive stress at which cro ss— se ct i onal distortion 
begins in columns with the sections considered in this pa- 
per. If, t and "b are the thickness and the half-width" 
of the flange, respectively i for an I-section, then the 
restraining effect of the web, whether positive, negative, 
or zero, is included in the coefficient k. On the other 
hand, if t and b are the thickness and the width of 
the web, a different set of values for k is obtained. 
It is therefore necessary to decide whether t and b in 
the equation for the critical stress shall refer to the 
flange plate or the web plate. Each form has advantages 
in certain cases. In this report^ both forms" will be 
given, either of which may be used to calculate the critical 
st res s . 

As applied to -I-section columns, these two forms are: 

f = 1 _ - - ( 2 ) 

cr 12(1 - n 2 )b3P a .- 

and 

12(1 - (X )b w 

In these formulas, either of which may be used, 

tjji and t^y are the thicknesses of the flange and_the_ 

web plates, respectively. 
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ky and k^, nondimensi onal coefficients dependent on 

the shape of the cross section. (See 
f igs . 3 and 4. ) 

"Dp is the half -width of the flange. (See 
sketches on figs. 3 and 4.) 

b^, width of the web. 

The values of ky and k ff , as calculated "by the previ- 
ously mentioned energy method, are list-ed in tables I and 
II. These calculated values were plotted against h^/by 

and hp/h^, respectively, to obtain the curves shown in 
figures 3 and 4. 

Equations (2) and (3) are given on figures 3 and 4, 
respectively, with f cr replaced "by f cr /T). The nondi- 

monsional coefficient T) allows for the effect of stress 
beyond the elastic range and is discussed in a later sec- 
tion. 

-The relation betTre"e~n k-g. and k^ for a given I- 

section column is sotaetiines of interest; it is obtained 
"by equating the right-hand sides of aquations (2) and ( 3 ) , 
which gives 



\bTy/> ^-tji 

Equations (2), ( 3 ) , and (4) also apply to columns of 
channel- secti on and Z-section, for which the values of 
kp and k^ are given in figures 5 and 6. It is impor- 
tant to note that, in these cases, by refers to the to- 
tal width of the flange. The values of kp and kw used 
to draw the curves in figures 5 and 6 are listed in tables 
III and. 17, respectively. 

For the case of the rectangular tube, let the cross- 
sectional dimensions of the tube be as shown in the sketch 
on figure ?. It will always be possible to- select the 

sides of width b and h so that ^^1- With this 

choice of sides, the critical compressive stre-ss for the 
tube is 
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kTT G St h a 



(5) 



f 



3 



12(1 - \i )h 



in which k may "be read from figure 7. The values of k 
used to draw these curves are listed in table V. 

The k-coef f icients ^iven herein apply to columns in 
which the material is "both elastic and isotropic. Steel, 
aluminum alloys, and other metallic materials usually sat 
isfy these conditions provided that the material is. not 
stressed "beyond the elastic range. TThen a material is 
stressed "beyond the proportional limit in one direction, 
it is no longer elastic and is probably ho longer isotrop 
ic. In a later part of the papor a method is presented 
for calculating the critical stress when the columns are 
loaded beyond the proportional limit. 



The plate elements that make. up the sid.es of the col 
umns treated fall into two classes, namely; 

(a) Plates restrained along both edges, as the web 
of a channel or as the sides of the rectangular tube. 

(b) Plates restrained along only one edge, as one 
flange of the channel section or the half-flange of the 
I-section. 



The deflection equation assumed for plates of class 
(a) is " 



DEFLECTION EQUATIONS 




and, for class (b), 





(7) 



where 
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w is deflection normal to plate. 
A, B, 0, and D, arbitrary deflection amplitudes, 
y, coordinate- across plate, measured from one odgo. 
b, width of member. 

n, -number of half-waves in length L. 
x, coordinate in direction of length L. 
L, length of member. 

For the channel* the Z-, and the I-sections, the val- 
ues of 0 and D for the flanges may be expressed in 
terms of A and B for the web through the conditions 
that the corner angles are maintained during buckling and 
that the moments at each corner are in equilibrium . The 
values of B/A and l/n are thon adjusted to make the 
critical stress a minimum. The same considerations apply 
to the rectangular tube, except that two eq.uati.ons of the 
form of (5) are used, one for each pair of walls. 

The foregoing deflection equations used in the- energy 
solution were carefully selected. A comparison of the ox- 
act values of k with the values of k for the rectangu- 
lar tube as gLven by the energy method reveals that the en- 
ergy values are less than 1 percent in error. (See refer- 
ence 3.) The values of k for the other sections are also 
"believed to tie correct to the same order of magnitude. . 
This belief is justified because, in the limiting cases 
for which exact solutions are available, the precision is 
within these limits. In addition^ other problems in which 
these deflection equations have been used gave a precision 
better than 1 percent. 

DISCUSSION OF OHAHTS 



Figures 3 to 7 give the computed values of the k-coof- 
ficients plotted against b^/b-g,, bj/b^, . or b/h. When 

the webs are narrow in comparison with the flanges, the in- 
stability occurs first in the flanges. As b^/bp In- 
creases, a point is reached where the webs be.como the 
weaker part of the cross section. For example, this point 
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is clearly di scernible in figure 3 where two of the curves 
"break sharply at certain values of "b^/bj,. 

CRITICAL STRESS FOR LOADING BEYOND THE PROPORTIONAL LIMIT 



In the elastic range, the critical compressive stress 
for an ordinary column that fails by "bending is .given by 
the Euler formula. Beyond the proportional limit marking 
the upper end of the elastic range, the reduced slope of 
.the stress-strain curve requires that an effective modulus 
E be substituted for Young's modulus E in the Euler' 
formula. The value of E is sometimes written 

E = TE (8) 



The value of the nondimensional coefficient T varies 
with stress. By the use of the double-modulus theory of 
column action, theoretical values of T can be obtained 
from the compressive stress-strain curve for the material 
(reference 6, p. 572, and references 7 and 8)_. Tests show 
that, in practice, theoretical values of T, derived, on 
the assumption that no deflection occurs until the criti- 
cal load is reached, are too lar§e. It is theref o re^Jbe st , 
for practical use, to obtain the value of T from the ac- 
cepted column curve for the material, as will be shown" 
later. The values of T thus obtained take into, account 
the effect of imperfections that cause deflection from the 
beginning of loading as well as other factors that may 
have a bearing on the strength. 

Eor cross-sectional distortion of a thin-wall column, 
the critical compressive stress is given by the basic equa- 
tion (l). Beyond the proportional limit, the critical com- 
pressive stress is given by this equation with an effec- 
tive modulus TIE substituted for E, or: 

f cr = Tl -^^4~s .. ....(9) . 

12(1 - \i )b 

In the absence of. adequate test data, the_ value of the non- 
dimensional coefficient T| cannot be definitely estab- 
lished. It is reasonable to suppose, however, that Tl 
and T are related in some way. 



"Various equations relating T| and T have been sug- 
gested. The discussion of reference 3 points out that, 
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when T) Is considered to he a function of T, the equa- 
tion for T) \7ill depend upon' the manner of - evaluati on of 
T. If T is determined from the st re s s- st rain curve on 
the assumption that no deflection takes place until tho 
critical load Is reached, the effect of deflections from 
the beginning of loading must be separately considorod. 
If T is determined, however, from the accept-ed column 
curve for the material in the manner outlined in the illus- 
trative problems of reference 3, part, if not all, of this 
effect is automatically considered. 

A careful study of- the theory and of such. experimen- 
tal data as are available indicates that a conservative 
assumption is 

n = I-tr|j^ (10 ) 

provided that T is evaluated by use of the accepted col- 
umn curve for the material. This equation will probably 
need to be modified as more test data become available. 

Now t is iirself a function of the stress f cr . 
Hence, T) is a function of f cr * Consequently, equation 
(9) cannot be solved directly for f cr . If the equation 
is divided by T) , however, then f cr /T| is given directly 
by the geometrical dimensions of tho cross section. 

Thus, for the I-section column, equations (2) and 
(3) become 



fcr fc^Et/ 



and 



12(1--- V? )b ff 2 



£cr = __V^*w" 

1 12(1 - v. )b w 2 



(11) 



(12) 



For the columns of channel and Z-section, equations (ll) 
and (12) also hold, except that bji is thro entire width 
of the flange. 



For the rectangular tube 
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£g. - fc "-'\ . ■ - (is) 

1 12(1 - n )h 

The values of the k-coeff icionts are read as usual 
from figures 3 to 7. Thus, f CT /^ is determined for the 

cross soction under consideration from equation (ll), (12) % 
or (13). The value of f c is' then determined from 

RELATION BETWEEN f cr AND f aT /f\ 



Assume that the value of Tj is given oy equation (10 ), 
As stated previously, T and hence T) depend upon the 
critical stress. Although theoretically the values of f 
and hence of Tl can "be obtained from the stress-strain 
curve, they are "best obtained from the accepted column 
curve for the material. 

The equations that show the variation of T with 
stress for 24ST aluminum alloy that just meats the require- 
ments of Navy Department Specification 46A9a" (tensile yield 
strength, 42,000 pounds per square inch) are given in part 
I of reference 9. In order to show how similar equations 
can he derived for any other material, equations will he 
derived from the column formulas- given in reference 9 for 
24ST aluminum alloy. 

The accepted column formulas for 24ST aluminum alloy 
are given by equations (8) and (9) of reference 9. These 
equations are: 

For. 41,200 > f cr > 19,600 lb./sq. in.., 

f cr = 43,700 (l - Q. 00752 (14 > 

For f cr < 19,600 lb./sq. in., 

* _ 105200000 ." _ ... /-,*\ 

For the same member, the critical stress given by the 
accepted column curve .must equal the critical stress given 
by the Euler formula with an effective modulus E = TE 
substituted for Young's modulus E, or 
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- tt2tE . j - T 105200000 (16) 
cr " (L/p7 " (L/P) 3 

If equations (14) and (15) are solved for L/ p , the 
following expressions are obtained: 

For 41,200 > f cr > 19,600 lb./sq. in., 

. 43700 - f cr 
L /P - .28.6 — ' 

For f cr < 19,600 lb./sq. in. 



J-/ 



105200000 



(IS) 



Substitution of these values of L/p in equation (16) 
gi ve s ; 

For 41,200 > f cr > 19,600 lb./sq. in., 

fx fx 

1 cr A oo„ . _er \ 



T - 



A >?2 4 - (19) 
V 25700/ 



8925 

For f cr < 19,600 lb./sq. in., 

T = 1 (20) 

By the use of equations (19) and (20), the value of 
T can be established for assumed values of f cr « T ^- Q 
values of T obtained can then be substituted in equation 
(10) to obtain the corresponding values of T). If the as- 
sumed values of f cr are divided by the corresponding 

values of T|, a curve of f cr against Z CT /f\ can be 
plotted. In figure 8, several such curves are given for 
24ST aluminum alloy for different assumed relations be- 
tween T] and T. When the valuo of f cr /T| has beon ob- 
tained by use of equation (ll), (12), or (13), the value 
of f cr is read from, one of the curves of figure 8. 



The ultimate strength of a thin-wall column will-, in 
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general, be higher than the load at which cross-sectional 
distortion "begins. At stresses approaching the yield 
point of the material, the critical load and the ultimate 
load approach the same value. No attempt has he en made in 
this paper to discuss the ultimate strength of a thin-wall 
column; the solution for the critical load logically pre- 
cedes the solution for the ultimate load. 



ILLUSTRATIVE PROBLEMS 



It is desired to calculate the critical compressive 
stress at which cross-sectional distortion "begins in two 
I-section columns constructed of 24ST aluminum alloy: 

Column A C olumn B 

"bp = 2 in. 

"b w = 4 in. 

t F = 0.2 in. 

t^ = 0.1 in. 

Solution for Column A 



^1 - - - 2 

Vff ~ 1 

tj. 0.1 

k ff = 0.662 

3 = 10.66 x 10 6 lb. per sq.. in. 



bji = 1 in. 

= 2 in. 

t-g, = 0.1 in. 

t« = 0. 1. in. 



Prom figure 3 , 
For this material, 



and (j, = 0.3 

Prom equation (ll), 



% = 0.662 X ^ X 10. 66^10* x .(O.l) 3 . Q 
n 12(1 - 0.2T) X (l) 3 

Prom the 'solid curve of figure 8, f = 33,000 lb./sq.. in. 
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Solution for Column B 

-I = - = 0.5 
b W 4 

*w = 0.1 = 0>5 

tj O.R 
From figure 4, k^ = 6.82 

From this material, E = 10.66 x 10 6 Id. per eq.- in. 
and \i = 0.3 

From equation (12), 

f e _ 6.82 X tt 3 X 10.66 x 10 S X (0.1) 8 y 

-~ = r -ST -Ta = 41,100 Id./scl. in. 

n 12(1 - 0.3 a ) x (4) a 

From the solid curve of figure 8, f cr = 28,800 lb./sq. * n ' 

CONCLUSIONS 



1. The critical compressive stress at which cross- 
sectional distortion begins in a thin-wall column of- chan- 
nel, Z-, or I-section is given by either of the following 
equat i ons : 

k^Et-p 8 

= Tj 



or 



cr ~ 'i : a. 2 

12(1.- u- )b F 



f m V^V 



cr ~ a. a 

12(1 - U- )bw 

where E and (J- are Young's modulus and Poisson's ratio 

for the material, respectively. 



bj, is the half-width of the flange for the 
I-section, or the total width of the 
flange f_or the channel and Z-sections. 
(See sketches on figs. 3 to 6.) 
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b^, width of the web. 

tj, and t w are the thicknesses of the flange and 

the web plates, respectively. 

kji and k w , nondimensional coefficients read from 

figures 3 to 6. 

T] is a nondimensional factor taken so that 
TIE gives the effective modulus of 
flange and web at stresses "beyond the 
elastic range. 

For the rectangular 'tube , 

k-rr^t^ 3 - - - - 

f cr = ^ — : s—~s 



where 

h 



12(1 - v )h 



and t^ are the width and thickness, respective- 
ly, of the long side of the rectangu- 
lar-tube section. 

k, nondimensional coefficient road from 
figure 7. 



2. At stresses "beyond the elastic range, the value 
of the effective modulus TIE for local "buckling "of thin- 
wall columns will depend upon tests. In the absence of 
such tests, however, it is reasonable to assume that Ti 
is a function of T, where TE is the effective modulus' 
of an .ordinary column at stresses beyond the elastic range, 
A careful study of theory .and of such experimental data as 
are available indicates that it is safe to assume 

T + 3 

i ; 

provided that T is evaluated by use of the accepted col- 
umn curve for the material. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Ta. , September 13, 1939. 
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TABLE I 

Calculated Minimum Values of 
kj, for I-Sections "by the Energy -Solution 



(H = 0.3) 



X. tw 


k E 


X tin 


0.5 


1 


2 


X 

.^k 








U 


1.288 


1.288 


1.288 


. U 




1 . 019 




. 4U 


.623 




1.193 


a D U 




. 852 




• HU 


.567 




1.147 






. 783 




. 1.20 


.5 36 




1.127 • 


J. • ffcU 




.735 




1.60 


.508 




1.118" 


1.70 


.500 






1.80 


.491 


. 688 




1.85 


.486 






1.90 


.473 






2.00 


.427 


.662 


1.113 


2.40 


.295 


.597 


1.110 


2.80 


.216 


.520 


1.105 


3.20 


.165 


.440 


1.093 


3.40 






1.062 


3.60 


.130 


.367 


1.015 


3.80 






.955 


4,00 


.105 


.306 


.893 


5.00 


.067 


.200 


,623 


6.67 

- 


.037 


.112 


.364 
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TABLE II 
Calculated Minimum Values of 



for I-Sections "by the Energy Solution 
(jx = 0.3) 











"b w \ 


0.5 


1 


2 


0 

,050 
.100 ■ 
.150 
. 200 
.250 
.278 
.313 
.357 • 
.417 
.5 00 * 
.526 
.556 
.588 
. 625 
.714 
.833 
1 . 000 -. 


4.00 
6.05 
6.46 
6.61 
6 6 8 
6.72 
6. 73 
6.75 
6.77 
6.79 
6.82 
6 . 84 
6.37 
5.78 
5.20 

3.09 


4.00 
4.49 
4.82 
4.98 
5 01 
4.89 
4.75 
4.50 
4.08 
3.44 
2.65 

2.23 
1.44 ■ 

— • 

.784 


4.00 
4.06 
4.09 
4.04 
3.89 
3.57 
3.29 
2.80 
2.17 
1.60 
1.11 

.715 
.406 
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A 

TABLE III 



Calculated Minimum Values of . 
for Channel and Z-Sections "by the Energy Solution 



LiL. = 0-5) 





k E 


"b-gt N. t-ci 
H N. £ 


0.5 






1 


2 


°r \ 






(J 


1 P R R 
J. I coo 


1 PAR 
J- « c a a 


1 P 8R 

1 i cuo 


. 200 




T ITT 
X . ill 




.400 


sac 
• o y o 




X » di «a*fc. . 


.6 00 




Q A P 
■ Sue 




. 800 


APT 
, Dei 






1.00 




RQ P 
• o 3 e: 




1.20 . 


lUfO < 




J-i 1 JO 


1 . 40 




R ^ A 




1 . oO 


■ uco 




X • X O O 


1.75 


*5 OA 






1,80 


.499 


.770 


- 


1 DO 


. 49 3 






1.90 


.455 






p on 


.410 


. 730 


1.187 


2. 20 


.338 






2.40 


.284 


. 529 


1.188 


2.80 


.208. 


.521 


1.190 


3.20 .. 


.159 


.423 


1.192 


3.40 






1.178 


3.60 


.125 


.345 


1.103 


3.80 






1.021 


4. 00 


.101 


. 284 


. 340 


4.40 


.083. 


.236 


.799 


4. 80 




.199 


.681 


5.20 


.059 


.170 


.587 


5 .60 




.146 


.508' 


6. 00 


. 044 


.127 


.444 
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TABLE IV 

Calculated Minimum Values of k w 
for Channel and Z-Sections "by the Energy Solution 



(H = 0.3) 



\ 4* 

11 




W 






0.5 


1 


c. 


U W \ 






- 


0 


4.00 


4.00 


4.00 


. 050 


5.46 


4.26 


4.03 


.100 


6. 02- 


4.45 


4.04 


.130 


6.19 


_ 


_ 


.167 


6.31 


4.58 
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TABLE V 

Calculated Minimum Values of k 
for Rectangular Tube "by the Energy Solution 
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